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Abstract 

• i-H , 

' We investigate the role of a certain class of QCD corrections to electroweak vector 

^ , boson scattering at small scattering angles and large energies. These are present since, 

from the perturbative analysis, the vector bosons may dissociate into quark- antiquark 
pairs giving rise to colour dipoles interacting through gluon exchanges. After the 
computation of the vector boson impact factors, we present expressions for the lowest 
order QCD scattering amplitude and for the leading logarithmic BFKL amplitude. 
Particularly we discuss numerical results for the process 77 — ► ZZ. The QCD correc- 
tions to the cross section resulting from the interference with the electroweak ones are 
estimated and compared with the leading pure electroweak part. Corrections resulting 
from the leading log BFKL amplitude are of the order of few percent already at the 
0.5 — 1 TeV energy range. 
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1 Introduction 



The search for physics beyond the Standard Model (SM) and the determination of the 
crucial and still not understood mechanism of the electroweak symmetry breaking in the 
SM itself are major subjects of investigation in particle physics. There is the hope to in- 
crease the understanding with the help of the next generation of supercolliders, and major 
efforts in this direction have currently been made. First insights are expected from the 
future large hadron collider (LHC) runs. On the other hand high precision measurements 
at the Next Linear Collider (NLC) are considered to be of great importance, with features 
complementary to LHC. Next Linear Collider is a generic name for a e + e~ machine oper- 
ating in the energy regime up to 1 TeV, providing a very clean environment, and making 
available observables for many processes at an accuracy better than the percent level. In 
addition to the e + e~ collider mode, one may have the capability of running the NLC in 
a 77 collision mode via Compton backscattering of laser photons off the linear collider 
electrons. 

One of the important processes one may consider in a linear collider is the scattering 
of vector bosons. Indeed, at the NLC the electron and positron bunches can radiate vector 
bosons which interact thereafter. Diverse scenarios are considered also in the Compton 
collider mode. Each incoming photon may turn into a virtual WW pair, followed by the 
scattering of one W of each pair to end with WW or ZZ states pp. There are also the 
processes where initial photons scatter directly into arbitrary vector bosons. 

Vector boson scattering is related to the open question of mass generation. Since 
the longitudinal modes of the massive vector bosons reflect the nature of the Goldstone 
modes of the unbroken theory, these may give important informations about the nature 
of electroweak symmetry breaking. Without the discovery of a scalar Higgs particle the 
vector bosons have to interact strongly at TeV scales, in order to preserve unitarity. In 
this scenario the vector boson scattering will be one of the key processes to be studied in 
the nearest future j2J. 

Another important question one may pose while searching for new physics is the ex- 
istence of anomalous triple and quartic vector boson couplings which may affect vector 
boson scattering. Contact interactions are regarded as a possible extension of the SM, 
which can be well investigated with vector boson scattering . The natural order of mag- 
nitude of the anomalous couplings [I] is small and, in order to eventually extract them, one 
needs to both measure and compute theoretically the SM cross sections with a precision 
better than 1%. 

In this context the 77 — ► ZZ process has a special role. Since there is no perturbative 
contribution at the tree level, such a process is sensitive to the particles and the new 
physics phenomena which contribute through radiative corrections 0- This leads to an 
independent and complementary method to the direct production of new particles. In 
the same process the detection of CP violating phases |H| or additional effects due to 
the exchange of Kaluza-Klein gravitons in large extra-dimension scenarios have been 
discussed. 

In summary, vector boson scattering processes may be seen as a central tool to probe 
physics beyond the SM and to reveal the nature of symmetry breaking. In order to exploit 
these features, a high precision is needed, both on experimental and on theoretical side. 
Calculations of these processes at the lowest available order may not be accurate enough. 
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One loop corrections have been calculated for all those vector bosonic processes [S] where 
a tree level contribution is present. For the 77 — > ZZ process, which starts at one loop 
in (Xewi 110 higher order corrections are available yet. The question naturally arise, if this 
accuracy is high enough to disentangle new physics from the SM, providing a window to 
non-accessible energy regimes. 

The aim of this paper is to address the question if the QCD corrections could play a 
role in vector boson scattering processes. Naively this seems not to be the case, since they 
arise only at higher loop corrections. On the other hand, even if the lowest order non- 
suppressed QCD correction to the amplitude for this process is smaller by 0(a 2 s ) w.r.t. the 
electroweak one, in the cross section the interference among the two can reach the percent 
level. Moreover at high energies there are kinematical regimes wherein the emission of 
soft gluons generates large logs in the energy which make the QCD amplitude rising with 
energy. 

The electroweak cross section of VV — ► VV (with V = W,Z, 7) gains its biggest con- 
tribution from small scattering angles, and with rising scattering energy the cross sections 
are approaching a constant in computations up to the one loop level. This contribution 
is coming from diagrams where only spin-1 particles are exchanged in the t-channel. At 
leading order and 0(a 2 EW ) these exchanged particles are only electroweak vector bosons. 

The same configuration arises for QCD corrections with gluons in the t-channel. The 
first gluonic contributions, which arise at the two loop level, are suppressed since these 
gluons are always accompanied by fermions in the t-channel. The non-suppressed QCD 
contributions occur at higher orders. These emerge when the vector bosons fluctuate into a 
quark-antiquark pair and these dipoles interact through gluons. At the lowest order, when 
two gluons are exchanged in the t-channel, one obtains a contribution constant in energy 
while higher order corrections provide large logarithms in the energy. These contributions 
cannot be neglected at large energies. The extreme attitude is to consider a resummation 
described by the BFKL equation [2j which gives an upper bound estimate of these effects. 

We study in detail these QCD corrections for the reaction 77 — ► ZZ, whose cross 
section up to now was only studied at leading order (one loop in the electroweak sector) 
by G. Jikia ^Oj. His results show that the cross section is approaching a constant for 
high energies and its main contribution comes from the helicity conserving amplitudes at 
small scattering angles. This is exactly the kinematical region where the above discussed 
QCD contributions have their biggest influence. 

We begin with the analytic calculations for an arbitrary bosonic process. In the next 
section the two gluon exchange amplitude is computed, including the calculation of the 
non- forward impact factor. Since we consider the high energy limit the amplitude factorise 
into two independent impact factors, built from the bosonic wave functions which describe 
the probability that a vector Boson fluctuates into a quark anti-quark pair. In the third 
section the cross section of the resummed gluonic contribution is presented, i.e. the BFKL 
amplitude. In order to give some numerical estimates, we concentrate on the definite 
process 77 — ► ZZ . In the fourth section we calculate the EW and QCD parts of the 
amplitude in the Regge limit. Next we evaluate the corrections numerically. Because of the 
strong scale dependence we can give only estimates of the corrections, which we find to be 
of the order of a few percent. Finally we present our conclusions. Our result suggest that, 
in particular in the small angle region, possible deviations from EW calculations cannot 
be interpreted as signals of new physics without taking into account the QCD corrections. 
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Figure 1: The two gluon exchange of the VV — > W elastic cross section. 



2 Impact factors and lowest order QCD contribution 

At high energies the lowest order leading QCD contribution to the vector boson scat- 
tering is pictured in FigQ where the lowest order diagrams with two-gluon exchange are 
illustrated. The vector boson fluctuation into a quark anti-quark pair is described by the 
product of the two boson wave functions. This quantity is called the (non-forward) vec- 
tor boson impact factor. Due to the high-energy factorization one can compute first the 
impact factors and then integrate these with the gluon propagators. In the high-energy 
limit, corresponding to s going to infinity and t fixed (Regge limit), we calculate first the 
s-channel discontinuity and restore after the amplitude through dispersion relations. Thus 
the intermediate quark lines inside the impact factors are taken on-shell. This procedure 
and the high-energy factorization simplifies the calculation significantly. The s-channel 
discontinuity of the scattering amplitude takes the form: 



where <3? is the impact factor and the boldface marked letters are the transverse part of 
the euclidean momenta. As usual, the integration of the two longitudinal components of 
the loop momenta are absorbed into the definition of the impact factors. The remaining 
integrals may be computed numerically. 

The momenta of the incoming bosons are q + I and p — I, while the outgoing bosons 
have momenta q — I and p + l. We use Sudakov variables with the light cone vectors 
p' = p + yq', y ~ —p 2 /2p' ■ q' and q' = q + xp', x ~ —q 2 /2p' ■ q', while s = (p + q) 2 ~ 2p' ■ q' 
denotes as usual the squared center of mass energy of the VV scattering process. Then 
one parametrizes any 4-momentum u as 



The impact factor for incoming and outgoing photons has been calculated recently 
and since we adopt the momentum assignment which was given there, we refer the reader 
to this paper for further details of the calculation. 





u = a u q' + (3 u p' + u± . 



(2) 
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The incoming and outgoing bosons are characterized by the invariants Q\ = —(q + 
I) 2 , Q- = —{q — 1) 2 and by similar expressions associated to the other two bosons involved 
in the process. For the full process the Mandelstam variable t = (2l) 2 has to be considered. 
In the case of on-shell bosons the virtuality reduces to the boson mass, i.e. Q 2 = —My. 

We remind that in the limit of Regge kinematics one can neglect terms which are sup- 
pressed by powers of t/s and Q±/s, leading to important simplifications. Let us consider 
the polarization vectors for the upper incoming and outgoing bosons. Their longitudinal 
polarization in terms of the Sudakov variables are given by 



e L (q±l) 



1 



IQ 



± 



/ 2I 2 \ 

[l± ai )q'+ zTA-y \P'±1- 



while the transverse polarisation vectors are 



(g±0 



Ah) 
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Ah) 



(3) 



(q' -p'±l±), 



(4) 



where h = ± denotes two helicity states and 



Ah) 



-^(0,l,±i,0). 



(5) 



Again, for the bosons appearing in the lower impact factor the same expression holds with 
an exchange of q and p and I —I. 

A difference between the photon and the other bosons is given by the chiral coupling 
to the quarks of the latters. Consequently extra terms in the traces appear and the 
ones proportional to odd powers of 7 5 do not contribute because of charge conjugation 
invariance. Thus, regarding the impact factors, the W and Z bosons behave both in the 
same way as the photon, i.e. we can use the result from and only the magnitude of 
the coupling to the quarks depends on the nature of the boson. 

The impact factor where both of the bosons have longitudinal polarisations takes the 
form: 



$LL(q,k,l) = a s ^N 2 - 1 ^ C / 4v / 4^ \Q+\\Q-\ J da J 



d 2 u 



xa 2 (l — a) 



while for transverse polarized bosons we have: 




Do 



(6) 




Do 



(7) 
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Finally, if the incoming boson is transverse polarized while the outgoing one is longi- 
tudinal, the impact factor reads: 



*gj.( S , k, () = a./NZ - 1 YJ C, 2Vfe / da / ^ x (8) 




X a (l -*)(!- 2.) I^-J+J f^-^j-«f, 



and viceversa for the incoming boson longitudinal and the outgoing one transverse, one 
has 

fc, I) = a 8 yjN* - 1 ^ Of 2v^ |Q_| ^ da | (9) 
x a(l — a)(l — 2a) ei • 















* 


^2) 



where e, = 1/v 2(1, ±z) are two-dimensional polarization vectors corresponding to the two 
transverse polarizations j = ±. We have also used: 

Df = (u± (1 - a)l) 2 + a(l - a)Q 2 ± + mj 

D± = (u - k =F al) 2 + a(l - a) Q± + m| 

= u±(l-a)l 

= u-k=Fal- (10) 

Marked in boldface are the Euclidean form of the transverse momenta, i.e. u\ = — u 2 < 0. 
The constant Cf is the product of the coupling of the incoming and outgoing boson to 
the quarks of flavor /. The impact factors have been written in a factorised form where 
the two parts correspond to the boson wave functions, describing the probability that a 
boson dissociates into a quark-antiquark pair. In order to obtain the lowest order QCD 
amplitude, one has to substitute these impact factors into eq.([Tj) and integrate over the 
gluon momenta. 

The two helicity-flip impact factors are in general non-zero. They vanish only for 
forward scattering, where 1 = 0. This is due to the different symmetry behaviour of 
the transverse and longitudinal boson wave functions. This property is more trasparent 
in the coordinate space formulation, eq.(34) of [TTj. Here the impact factor is written 
as a convolution of the two boson wave functions and the dipol interaction. The dipole 
interaction is invariant under the transformation of the dipole size vector r — > — r, and 
so is the longitudinal wave function, eq.(22) of |llj . In contrast to this the transverse 
wave function, eq. (28-31) of is antisymmetric under this transformation, since it is 
proportional to r. The convolution of two transverse wave functions is therefore symmetric 
under this transformation, whereas the convolution of a transverse and a longitudinal one 
is antisymmetric, leading to the vanishing result for 1 = 0. If 1 is non-zero, this symmetry 
properties are broken, resulting in a non- vanishing §tl and &tl- 
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Figure 2: The BFKL Pomeron exchange in the VV — > VV elastic cross section. 



3 Leading log (BFKL) contribution 

At high energies and close to the forward region (Regge limit), higher order contributions 
in a s start to be important. Diagrams of ladder topology, built with (non elementary) 
reggeized gluons, can contribute up to the order a™ln n (s/so)- On retaining the leading 
logarithmic (LL) contributions and performing a resummation a rough estimate (overesti- 
mate) of the total QCD corrections is obtained. Nonetheless one expects to find the QCD 
contribution to the process in the window constrained by the simple two gluon exchange 
on one side and the resummed LL BFKL Pomeron exchange on the other. 

This possibility is related to the property of the factorisation of the QCD scattering 
amplitude, in the Regge limit, in terms of the impact factors associated to the external 
particles and the BFKL Green's function, as illustrated in Fig. [21 

For the physical process we want to study the impact factors have been calculated 
in the preceeding section and the discontinuity of the LL resummed amplitude is simply 
obtained by inserting the non-forward BFKL Pomeron Green's function in the place of 
the two gluon propagators. Indeed one has 

Disc s A = is J -^^L$(g, k, l)G{y\k + l,k-l;k' + l,k' - l)$(p, —k', -I), (11) 

where y = In s / sq is the rapidity variable and sq is a typical scale in the process. Clearly, 
for a s — > 0, when all rungs of the BFKL resummation decouple, the Green's function 
reduces to describing the two gluon exchange, which means that we use the following 
normalisation: 

lim G(y\k + l,k - l;k' + l,k' - I) = <5 (2) (fc - k') (12) 

«.-o vy| ; (k + l) 2 (k-l) 2 K J K 1 

In the following we briefly give the ingredients needed to construct the BFKL Green's 
function. The eigenfunctions E h ^ of the BFKL kernel, which appears in the integral 
Bethe-Salpeter like equation resumming all the LL terms, are well known in coordinate 
space, wherein their form is dictated by conformal invariance. They are given by 



w^H^)\*k)*' (13> 
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with rij = ri - r,j etc., h = (1 + n)/2 + iv, h = (1 - n)/2 + iv (h* = 1 — h, h* = 1 - /i). 
Let us note that standard complex notation for two-dimensional vectors is used on the 
right-hand side. The conformal weights h, decomposed in n (integer conformal spin) and 
v (real), label the principal series of unitary representations of the Moebius group. Fourier 
transforming to momentum space (we use the Lipatov's convention which assigns a l/(2-7r) 2 
to any coordinate integration) one finds |15| : 

d 2 ri d 2 r 2 



E^ik.M) = J ^J^ E h,h(ri,r2)^ k ^ +k ^ = c(x(fe 1) fe 2 )+(-irX(fe 2 ,fe 1 ) 

(14) 

The coefficient C is given by 



(-i) n - - 
C = ±—^hh(l - h)(l - h)T(l - h)T(l - h). (15) 

(47Tj z 

The functions X in complex notation can be expressed in terms of hyper geometric func- 
tions: 

_ (|)"" 2 (|)" 2 F (l - H, 2 - Mi -f) F (l - S, 2 - S; 2; -| 

(16) 

This analytic form does not contain any term of the type 5 2 (k\) or 5 2 {k,2) which are present 
in the coordinate representation eq. (|13|) . The impact factor of a colorless external particle 
vanishes for zero momentum of any gluon exchanged and therefore the delta-function type 
terms do not contribute. 

The Pomeron intercept has the form op(0) = 1 + x(^> n )> where 

( /l + l n l /l + l n l \\ N c a s 
Xiy, n) = Xh = a s I 2-0(1) - ip{ — h iv) - ip{ — iv) J , a s = , (17) 

are the kernel eigenvalues corresponding to the eigenfunctions in equations (|13|) - <|14f) . We 
shall use the Pomeron Green's function in the momentum representation. Starting from 
the form in coordinate representation ^H] one can perform a Fourier transform and 
making use of the Casimir operator properties of the Mobius group, one finally obtains 
for the non-amputated Green's function: 

G 2 (y\k 1 ,k 2 ; kykv) = (2vr) 3 j d^{h) e y *» N h x (2vr) 2 E^fa, k 2 )El Ji (k v , fe 2 ,) , (18) 

where we use the measure in the conformal weight space J dfi(h) = J2 n I ^ v with the 
following normalization factor iV^: 

(^ 2 + 

h [u 2 + {n- l) 2 /4] [v 2 + {n + 1) 2 /4] ' 1 ' 

The (2-7r) 3 factor in front of the integral comes from the normalization fixed by ea. (ll2l) . 

Let us finally rewrite the expression of the amplitude discontinuity in ea. (|lljl by in- 
serting in it the relation 1)18(1 . On defining the conformal weight representation for the 
impact factors according to 

MQJ)= J ^Hq,k,l)E hJl (k + l,k-l), (20) 
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together with the primed & h (q,l) computed with the conjugate eigenstate, one has 



Disc 8 A = is(27rf / dfi(h) N h * h (g, l)<t>' h {p, I) . (21) 



We stress that due to high energy factorization the energy dependence is encoded in the 
LL BFKL Green's function and the only other dependence in the momenta of the incoming 
particles is through the invariant external masses on which the impact factors depend. In 
the last expression the dominant contribution for large y can be easily extracted with a 
saddle point analysis for n = around v = and therefore reduces to 



Disc s A = is(2vr) 5 16e^ -g^ $ ( g , l)& Q (p, I) , (22) 

where xo = 41n(2)a: s is the leading BFKL intercept and a = 14£(3)a s . 

The forward limit presents a well known cusp behavior and it must converge to the 
result obtained on using a forward BFKL Green's function 

G F (y\k,k>) = (2n) 3 J d»(h) e »*^g(fc)£g* \k') , (23) 

which is written in the spectral represention with the forward BFKL amputated eigen- 
functions given by E^(k) = l/^^k^k* h ~ 1 . 

The saddle point approximations can be carried on in a way similar to the non forward 
case. We shall define 

< ) ^ = /(^^' fc '°)^2 (fe) ' (24) 

'(F) 

as well as a ^ (q) for the scalar product with the conjugate Pomeron state and write 
the discontinuity in the forward case as 

Disc^A = ls (2n) 3 e y^^^4 F \q)% (F \p) . (25) 

Let us note that in the transition from the nonforward to forward physics the perturba- 
tive QCD analysis is increasingly affected by the long distance interactions (contributions 
for all values of the impact parameter), since gluons are massless. This implies that some 
care is required to treat the general non- forward expression in this limit. The saddle point 
approximation for n = around u = leading to eq.(|22|) is accurate only for I down to 
values of the order of one. Strictly speaking, performing the saddle point approximation 
of eq . 1)211) in the region of very small I one needs to shift the saddle point which becomes 
I dependent. In our numerical calculation we have used eq.© down to |21| 2 = 1/2 Mf, 
and then we have extrapolated to the forward value, eq. ([25j) . 

A more precise analysis, avoiding the use of the saddle point approximation is of 
course possible but it is not needed for the order of magnitude estimate which is pursued 
here. By the way few checks have been done with an exact Green's function computation, 
confirming such an expectation. 
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Figure 3: The surviving contributions to the leading electroweak amplitude in the Regge limit for 
real 7 and transverse polarised Z. 



Finally, in order to evaluate the elastic cross sections we need to obtain the amplitude 
from the computed discontinuity. Instead of using the full dispersion relation technique, 
a simple evaluation compatible with the approximation adopted is given by 



as provided by the signature factor in the partial wave amplitude analysis. 

4 The process 77 — > ZZ in the small angle limit 

We shall henceforth concentrate on a definite process, 77 — > ZZ, and compute for it a 
numerical estimate of the QCD contributions. We remind that we are interested in the 
kinematical region of small scattering angles, i.e. scattering close to the forward region, 
so that one can do the calculation in the high-energy approximation suited for the Regge 
limit: terms which are suppressed by powers of t/s can be neglected, with a significant 
simplification in the computation. We shall in fact use the results presented in section 2 
and 3. 

4.1 The leading electroweak contribution 

Let us first analyze the pure electroweak contributions to the process. This process starts 
at one loop and was calculated by Jikia ^01 ■ In order to calculate the QCD corrections to 
77 — > ZZ one has to know the leading electroweak contribution at the amplitude level, add 
to it the leading QCD amplitude and take the modulus squared. The main contribution 
to the corrections will come from the interference term. In the high energy limit where t is 
kept fixed most of the contributions are suppressed by powers of s. The surviving graphs 
are pictured in Fig. |21 Also for these only the helicity conserving amplitudes survive, with 
both of the bosons transverse polarised. One may confirm this by taking the limit s to 
infinity and keeping t fixed within the expressions of ^U] . 

In order to simplify the calculation we compute first the s-channel discontinuity, i.e. the 
s-channel internal W lines are taken on-shell. Here the s-channel discontinuity coincides 
with twice the imaginary part of the amplitude. With the same assignment of the momenta 
already used in the last sections the imaginary part of the box contribution reads: 



lmA^ zz = 114^2 m ~ k? ~ MlA 6[(p + kf - Myr] A^ WW A^ ZZ , (27) 




(26) 
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where a summation over the W boson helicities is implied. Performing the loop integration 
and retaining only the leading power of s we get: 



ImA 



8vr 



Box 



s a 



em o 



In 



t 2 



AtM 2 , 



-t + Jt 2 



4tM$y 



-t 



t 2 



4tM 2 Vi 



(28) 



where |21| 2 a 
that ln(— s) - 



—t was used. Equation (|28|) yields immediately the real part by just noting 
In s — m. Thus the reconstructed (leading) full box amplitude reads: 



— s 



In 



7T 



8vr 



a 



em o 



In 



AtM 2 , 



-t + Jt 2 



AtM 2 , 



t 2 



* tM w, 



(29) 



The amplitude of the crossed diagram is obtained from the box contribution (|29|) by 
replacing s by u ~ —s. Adding the box and crossed contributions, at the end one is left 
with the imaginary part of the box diagram: 



A++++ 



i Im A B ox 



(30) 



Thus in the Regge limit the leading order electroweak contribution is purely imaginary 
(s 3> —t), and the cross section is constant in s. Considering the limiting behavior of 
ea.(|30|) in the region s ^> — t and for the case t = one obtains the expressions 

presented in [TU] , 



4.2 QCD corrections 

The amplitude of this process in the high energy limit is illustrated in Fig. |^1 The in- 
gredients needed have been discussed in sections 2 and 3. To proceed with an explicit 
calculation we have to specify the couplings Cf, which depend on the external bosons, and 
the virtualities appearing in the impact factors given in equations (JJJ) and ©. The start- 
ing expressions are for massive quarks. For practical purposes only the top quark mass 
has to be taken into account because of the large Z mass. The integrals which correspond 
to the quark loop have to be performed. The last integration we compute is the scalar 
product of the impact factors with the BFKL states eq. ((2*0]) and (|2~I|l . 

We calculate this process for all external bosons on-shell. In such a case we have 
identically zero and Q 2 _ = —M\. The coupling constants Ct - the product of the 7-quark 
and Z-quark couplings - are given by: 



Cf = 



Atto 



em 



2 sin 9 W cos 6 W 



{(If 



■Tf -2qj sin 2 



.}. 



(31) 



First we consider the scattering of the photons into transverse Z bosons. To proceed 
we have to solve the remaining integrations for the impact factor <S>tt in eq.(J7J). In order 
to compute the integral over the two dimensional euclidean vector u we employ a standard 
Feynman parametrization and rewrite eq.([7J) as: 



a. 



N 2 



Att / da 



d 2 u 

WT 2 



dx x 



(32) 
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4a (1 — a)[ej • u e*j • u — x(l — x)e{ • p a e| • p a ] — (u 2 — x(l — x)p 2 + m 2 A{ei ■ e* 



[u 2 + x(l — x)p 2 



xm 



m 



212 



4a(l - a)[e r ue*-u- x(l - i)e» • p b e* • Pb] 



x(l - x)pg + m^)(ej • e*) 



[u 2 + x(l — x)p 2 — xm + m 2 ] 2 



Here we have defined p a = — k — (1 — 2a)l , Pb = 2al and m = a(l — a)(M| +ze). Strictly 
speaking, at this stage of the calculation m is taken negative and subsequently, after having 
performed the integrations, analytically continued to positive values obtaining the above 
definition. The first fraction corresponds to the two diagrams where one of the gluons is 
attached to the quark and the other to the antiquark. The second fraction corresponds 
to the two diagrams where the two gluons are attached to the same (anti)quark. Each of 
them would lead to a divergent integral but, or course, the sum in eq, (|32j) gives a finite 
result. 

On performing the u integration we end up with: 



a. 



N 2 



/ 



/ — f dx 

2-7T Jo 



4a(l — a) xil — x) x 



€j ■ p a £ ■ Pa 



6i • Pb e, • Pb 



x(l — x)p 2 — xm + m 2 x(l — x)p 2 



xm + m 



+ <k 



+ 




x(l — x)p 2 — xm + m 



x(l 



x)p\ 



xm + m 



x(l — x)p 2 — m 2 
x(l — x)p 2 — xm + m 2 x(l — x)p\ — xm + m 




(33) 



The computation of the x-integration is straightforward. Due to the lenght of the 
result we present here only the limit mj — > which preserves a significant simplification 
compared to the massive case. The corresponding helicity conserving result with massive 
fermions is presented in the appendix. After the integration in the Feynman parameter 
variable, one obtains in the zero quark mass limit for the separate helicity states: 



(— ) 

TT 

(+-) 
TT 



(q,k,l) 
(q,k,l) 



(-+) 
TT 



(q,k,l) 



a 



In 




p 2 — m 



I — 4a(l 



p x h p\ 



a) x 



2 
Pb 



m 



1 + T2 ln 



p\ 



m 



-m 



'TT 



'(q,k,l), 



(34) 

(35) 
(36) 



(37) 



with (p x ,p y 



p. For the helicity conserving parts in the massless quark limit the 



remaining a integration was also performed analytically and is presented in the appendix. 
For the result with massive quarks a numerical integration has been performed. 
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In the same way one can calculate the scattering of real photons into longitudinal 
Z-bosons and obtain (in the mt — > limit): 



^ L (q,k,l) =a s JNi-lY / C f 2V^M z I'p a(l-a)(l-2o) 

/ 



fi • Pb 

vl 



In 



Pb-m 



-m 



fi • Pa 

Pa 



In 



-m 



(38) 



In order to compute the two-gluon or the BFKL amplitude one has to substitute these 
expressions into eq.(^Q) or eqs. (|20L 124)) . 



5 Numerical results 

In this section we present numerical results on the QCD corrections for ZZ pair production 
in 77 collisions. Here we consider full circular polarization of the incoming photons. 
The relative QCD contributions are presented for transverse Z bosons only, because in 
the high-energy approximation the helicity changing parts of the electroweak amplitude 
are supressed by powers of s and thus are negligible. Nevertheless we shall display the 
helicity changing part of the BFKL amplitude as the pure QCD contribution alone. The 
parameters aw = o/s w , a = 1/128, mz = 91.2 GeV, m t = 174.3 GeV and a s (so) have 
been used throughout the numerical computations. We have been neglecting all quark 
masses apart from the top quark mass in the numerical calculation, since they are neglible 
compared to the Z boson mass. The inclusion of the top quark mass reduces the QCD 
amplitude by 25%. 

In terms of the amplitudes (|3()|21|2H)1 the cross section reads: 

^ = tLi IAew+Aqcd12 (39) 

where pt = |21| is the exchanged transverse momenta. The angular distribution of the 
cross section is obtained by the relation: 

dt = S -fi dcosO, (40) 

with (3 = \J\ — AMg/s. For small scattering angles in the high-energy approximation the 
two incoming and two outgoing particles do not 'communicate' with each other because 

of the huge separation in rapidity. As a result, the amplitude with the helicity -H > ++ 

is equal to the amplitude with helicity H ► H — etc. Since the helicity conserving am- 
plitudes completely dominate the cross section, after the summation of outgoing helicities 

the result is the same for all the possibilities of incoming helicities, i.e. -H ► TT is equal 

to H ► TT etc. In the following 77 — ► Z?Zt stands for all these possibilities as well as 

for unpolarised scattering. 

In Fig. 4 the angular distribution of the leading electroweak contribution for transverse 
Z pair production for polarised 77 collision at ^/s 77 = 1 TeV is given. The solid line is 
the result of the calculation in the high-energy approximation, presented in section 4.1. 
The dotted line is the result of the exact calculation, presented also in ^01 for ^/s 77 = 500 
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Figure 4: Angular distribution of the leading electroweak contribution for transverse Z pair pro- 
duction in polarised 77 collision at y/s = 1 TcV. The solid line is the result of the 
high-energy approximation, the dotted line is the result of the exact calculation. 



GeV. For cos 9 ~ 0.9 the result of the high-energy approximation has a deviation of ~ 
15% from the exact calculation, for smaller scattering angles the deviation is decreasing. 
With rising center of mass energy the accuracy of the Regge calculation is continously 
increasing. 

Let us consider the QCD corrections. In the Born analysis, from the naive counting of 
the couplings (see Fig. [I}, one expects in the interference term a suppression, with respect 
to the one loop EW case, due to the factor a? and an additional suppression of a factor 
~ 4 due to the difference in the EW couplings. This would lead to corrections of less 
than a percent. At high energies in the Regge limit we expect an enhancement due to the 
appearance of some large logarithms, which may be resummed in the BFKL scheme. We 
stress that the resummed leading log QCD corrections at the lower energies we consider are 
overshooting what is expected from the real contribution, nevertheless they provide a first 
estimate. We remind, as discussed previously, that the BFKL resummation is evaluated in 
the saddle point approximation. Fig. 5 displays the differential cross section with respect 
to the exchanged transverse momenta at a center of mass energy of ^/s 77 = 1 TeV. The 
solid line is the pure electroweak cross section, and the dashed line is the one with the QCD 
corrections included, eq.(j2U). The cross section has a slight rise approaching the forward 
region. This slight rise comes from the electroweak part since it completely dominates the 
amplitude. As in the following it will be visible, for the QCD part of the cross section the 
rise approaching the forward region is much stronger. Therefore the QCD corrections are 
enhanced close to the forward region. 

This behavior is shown in Fig. 6, where the QCD corrections to the differential cross 
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Figure 5: The pure EW differential cross section (solid line) versus the exchanged transverse mo- 
menta at y/s = 1 TeV and the one which include the QCD corrections (dashed line). For 
the BFKL contribution the scale s = (l/2Mz) 2 was used. 

section relative to the pure electroweak contribution at center of mass energies of 500 GeV 
(a), 800 GeV (b), 1.5 TeV (c) and 3 TeV (d) are shown as functions of j4/M|. At these 
energies = 4 corresponds to values of cos 6 (where 9 is the scattering angle) of 0.75 

(a), 0.90 (b), 0.97 (c) and 0.99 (d) respectively. This relative correction is defined as 

a _ | d<TQCD+EW da EW \ doEW 

\ dp\ dp 2 J I dp\ ' 1 j 

Since for leading order BFKL the scale is not fixed, there is an uncertainty which affects 
the magnitude of the QCD corrections. The solid lines in Fig. 6 are results for the scale 
so = (1/2 Mz) 2 while the dashed lines correspond to sq = M§. At the scale sq = M| 
the minimal center of mass energy y/s, for which the saddle point approximation for the 
leading logarithmic calculations gives a reliable result, is 800 GeV. Corrections resulting 
from leading order BFKL contribution are at the percent level. For higher pt they are 
about one percent and approaching the forward region (smaller pt values) they rise and 
lead to corrections of a few percent. Close to the forward region the sign of the corrections 
is changing as a result of the fact that for small 1 the imaginary part of (eq.20) gets 
an enhancement. Fig. 6 (a-d) summarizes the dependencies of the corrections resulting 
from leading order BFKL calculations with respect to the center of mass energy y/s, the 
exchanged transverse momenta pt and the scale sq. For rising center of mass energy the 
corrections increase. The solid lines (so = (1/2 Mz) 2 ) in Fig. 6(a) and (d) show that the 
corrections may rise in this regime by a factor of 2 if the energy increases from 500 GeV 
to 3 TeV. The scale sq implicates also a strong dependence. In fact close to the forward 



14 




Figure 6: QCD corrections to the differential cross section relative to the pure EW contribution 
for different scales sq and center of mass energies y/s. The relative correction is defined 

aS A = C d<T Q CD + EW _ dtTEW \ / d<TEW 

dp'^, dpij, 'I dp'^, 
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Figure 7: The QCD corrections to the integrated cross section relative to the pure EW contribution 
for different scales. The integration region was choosen as 1/4 M| < pr < 4M|. 

region the corrections also grow by a factor of 2 if the scale is decreasing from so = M§ 
to so = (1/2 Mz) 2 - The biggest dependence is coming from the exchanged momenta or 
equivalently from the scattering angle. The corrections have a strong rise with decreasing 
scattering angle or pt and change sign close to the forward region. 

In Fig. 7 relative corrections to the integrated cross section (ctew+qcd — &ew)/&ew 
are plotted as a function of the center of mass energy. These are also at the order of 
percent level. The corrections rise with increasing center of mass energy. The magnitude 
of the corrections is again scale dependent and the rise is stronger for smaller scales, as 
expected. Since the calculations are done in the Regge limit we integrated eq.jSHJ in the 
region of 1 /4 M§ < px < 4 M§ . An integration down to px = would change the result 
by less then 0.1 %, since the sign of the corrections changes very close to the forward 
region, (see Fig. 6). If one integrates over all scattering angles the corrections become 
less important, but since the main contribution of the cross section is coming from small 
scattering angles the result would not be affected strongly. 

In the following figures some properties of the pure QCD part of the amplitude are 
discussed. In Fig. 8 the absolute value of the QCD amplitudes is displayed for the helicity 
conserving case. The solid line is the BFKL (resummed) contribution while the dashed line 
is the Born (two-gluon) contribution. At y/s = 1 TeV the gluon resummation is enhancing 
the pure QCD Born amplitude by one order of magnitude. The same enhancement is 
visible in the QCD corrections to the electroweak cross section. Both amplitudes are 
rising with decreasing p?, whereas for higher energies the rise of the resummed amplitude 
will be stronger as the Born one. 
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0.5 1 1.5 2 2.5 3 3.5 4 



p^/M| 

Figure 8: The pure QCD part of the amplitude versus the exchanged transverse momenta. The 
solid line is the BFKL contribution and the dashed line is the QCD Born approximation 
(two-gluon amplitude). 

In Fig. 9 the absolute value of the two helicity changing parts of the amplitude at 
i/s = ITeV are plotted. The most noticable property is that the helicity changing parts are 
suppressed by more than two orders of magnitude with respect to the helicity conserving 

part. The -H ► cross section has, apart from its magnitude, properties similar to 

-H > ++, for example its slope as a function of pr- The -H ► LL part is also rising 

when pt is getting smaller but in the region of px = 1/2 Mz it has an extremum and is 
steeply going to zero for p? = 0. This behavior comes from eq.©, where the integrand is 
antisymmetric in a <-> (1 — a) for 1 = 0. 

6 Conclusions 

In this paper we have presented formulae for the leading QCD corrections to vector boson 
scattering at small scattering angles. These emerge since, like the photon, all vector bosons 
fluctuate into quark-antiquark pairs, which are strongly interacting colour dipoles. The 
nonforward boson impact factors describe the probability that the boson fluctuates into 
a quark-antiquark dipole. First the Born and after the leading logarithmic resummed 
(BFKL) contributions have been investigated for an arbitrary process of elastic scattering 
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0.5 1 1.5 2 2.5 3 3.5 4 

p2/M| 

Figure 9: The same as Fig. [5] for the helicity changing parts of the BFKL contribution, 
of vector bosons. 

In particular the process 77 — ► ZZ has been considered in detail. Both, the elec- 
troweak and the QCD parts of the amplitude have been calculated in the high-energy 
approximation. The results are displayed for the differential and integrated cross section, 
comparing QCD corrections relative to the one loop elektroweak calculation. We discussed 
some properties of the pure QCD part of the cross section as well, including also nonlead- 
ing contributions. The general result of the numerical calculation is that at s = 0(1 TeV) 
the QCD corrections affect the cross section of the process by a few percents. The correc- 
tions are smaller for larger pt values and rise with falling pr- Moreover they rise with the 
center of mass energy as expected. Since the leading log BFKL contribution ends up with 
a noticable scale dependence this results should not be regarded as a precision calculation. 
Therefore, we interpret this result as an indication that, in particular in the small angle 
region, possible deviations from EW calculations cannot be interpreted as signals of new 
physics without taking into account the QCD corrections. The same reasoning is valid for 
the integrated cross sections, since the biggest contribution to the cross section is coming 
from scattering at small angles. 
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7 Appendix 



In this appendix we present further formulas for the impact factor calculated in chapter 
4.2. 

In the case of massive fermions ea. (|34j) has the form: 



(++) 

TT 



(q,k,l) 



a) 2 ] x 



Pi 



vl 



- m + 2 
m + 2 TOj / 



log 



-m 



Pa 



A fi 



2 A,, 



log 



— m 
-m — 



Pa + A a 



log 



Pi 



A h 



-m 



\Pfe + A 6 



log 



A:r g 



-m 



-to - 
-m 



■p2 + A a 



-m 



Pi + ^b 



log 
- log 



-m+Pa 



-m + pi - A a 
-m + p 2 + A a 

-to +pl - A b 
-m + pl + A b 

A a 



-m + p 2 a + A a 

-m + pl - A b 
-m + pl + A b 



(42) 



while ea.(|38|) becomes: 



*$ L (q,k,l) 



iV c 2 - 1 V Ct -^=M Z C da (a(l - a)(l - 2a) 



gi • Pb / Pt_ 



m 



Ai, 



e i ■ Pa Pa~ m 



2pl 

with 



A n 



log 
log 



-m-p 2 h - A b 
-m-p\ + A b 

-m-p\- A a 
-m-p 2 a + A a 



log 
log 



-to + pi - A b 
-m + pl + A b 

-m + pl- A a 
-m + p 2 + A a 



log 
- log 





(43) 
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m 2 — m 
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A a 
A b 



(p 2 - m) 2 + 4m 2 f p 2 , 
(Pb ~ m ) 2 + ^ m }Pb > 



(44) 



and p 2 , pi and m was already defined in section 4.2. Taking the limit mj — ► one obtains 
eq.([34j) and eq.(|38j) respectively. This formulas were used for the top quark contribution, 
where the remaining a integration was done numerically. 

In the case of massless fermions for the helicity conserving part the a integration may 
be performed analytically, giving the result: 



<I> 



(++) 

TT 



(q,k,l) 



a s ^N 2 -l Cf {2 (xi -x\ + x 2 - x 2 ) + 

+ (2 - 3xi + 3x 2 - 2xfj log(l - xi) + (3xi - 3x 2 + 2xf) log(-xi) + 
+ (2 - 3x 2 + 3x 2 2 -2x1) log(l - x 2 ) + (3x 2 - 3x 2 2 + 2xfj log(-x 2 ) + 



+ 



3a 3 



-3<r + 3ad + 2(T J log(-) + 
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+ (a + d) (2a d- (2 a 2 + ad + 2d 2 ) log(l + ^ J, (45) 



with 



a = U 2 + Mz + ie 

b = -4(1 2 + kl) - Mf - ie 

c = (k + 1) 2 

d = -Mf - ie 



-b ± y/b 2 - 4ac 
*l,2 = ^ . (46) 
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